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The ground state of the Holstein molecular crystal model of a single electron in one dimension is 
studied using the Global-Local variational method, analyzing in particular the total energy, kinetic 
energy, phonon energy, and interaction energy over a broad region of the polaron parameter space. 
Through the application of objective criteria, a curve is identified that simply, accurately, and 
robustly locates the self-trapping transition separating small polaron and large polaron behavior. 
Particular attention is given to the kinetic energy, which is compared quantitatively with other 
perturbative and non-perturbative methods. 

PACS numbers: 71.38.+i, 71.15.-m, 71.35.Aa, 72.90.+y 



I. INTRODUCTION 



Despite the extensive effort and creativity that has 
been applied to the polaron problem for more than half a 
century, there remain significant aspects of polaron struc- 
ture and behavior that have defied satisfying explanation 
and quantitative description. Central among these is the 
so-called self-trapping transition, manifested as a strong 
feature in the intermediate-coupling regime that sepa- 
rates states that are predominantly "small-polaron-like" 
from states that are predominantly "large-polaron-like" . 

Though this transition is expected to be smooth at 
finite parameter values, and though perturbation theo- 
ries exist providing asymptotically accurate descriptions 
of polaron structure and properties on each side of this 
transition, the intermediate-coupling region in which the 
transition lies has proven very resistant to analysis. 

Although the physically-meaningful self-trapping tran- 
sition is not a singular feature of the polaron landscape 
at finite parameter values, the self-trapping transition 
is related to singular behavior in the adiabatic limit. 
This singular behavior in the far reaches of parameter 
space appears to influence both weak-coupling pertur- 
bation theory and strong coupling perturbation theory 
at finite parameter values in that weak-coupling pertur- 
bation theory breaks down as this feature is ap- 
proached from below, and strong-coupling perturbation 
theory p|-|l3|] breaks down as this feature is approached 
from above. 
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Much as a discrete curve (e.g., the unit circle) may 
describe a limited domain of convergence for a series rep- 
resentation of a function though the function itself may 
be well-behaved across most of that curve, there is mean- 
ing to the notion of a discrete line describing the location 
of the self-trapping transition though the physical phe- 
nomenon we study may be smooth throughout the finite 
parameter space. However, lacking absolute knowledge of 
the underlying self-trapping phenomenon, we must neces- 
sarily proceed more empirically and attempt to infer the 
location of a broadly-meaningful self-trapping line from 
limited observations of polaron properties. From such 
observations on numerous distinct quantities over large 
regions of the polaron parameter space, we conclude that 
the self-trapping line of the Holstcin model in one dimen- 
sion is accurately and robustly described by the simple 
relation 



g ST = 1 + y/J/hto , (1) 

in which J is the nearest-neighbor electronic hopping in- 
tegral, to is the Einstein frequency, and g is the dimen- 
sionless coupling strength. By this we do not mean that 
every physical property exhibits transition behavior on 
this curve at all parameter values, but that collectively 
transition behaviors in all measured properties are con- 
sistently and quantitatively related to this single curve. 
This curve could already be inferred roughly in Figure 4 
of Ref. |Q, but here we are led to this functional form 
by direct, high-precision numerical study. 

Describing the self-trapping transition with accuracy 
clearly requires methods that are non-perturbative. In re- 
cent years, several methods have been developed that are 
capable of describing the intermediate-coupling regime 
with high accuracy over non-trivial ranges of adiabatic- 
ity; these include variational techniques ^Q-Q, clus- 
ter diagonalization 11^,E3 E4[, density matrix renor- 



malization group (DMRG) |p5| 26 and quantum Monte 



Carlo simulations (QMC) [|27pp3|] 

Elsewhere |L8| , we have made direct quantitative com- 
parisons among a number of these and other indepen- 
dent methods, demonstrating broad quantitative agree- 
ment among the best. Here, we use the Global-Local 
variational method on periodic lattices of 32 sites, sup- 
ported by low orders of perturbation theory on infinite 
lattices. Details about the numerical method can be 
found in Refs. |@. 

We focus on the 1-D Holstcin Hamiltonian for a single 
electron 

H = Hkin + Hph + Hint , (2) 

Hkin = — J ^ a n {a n +i + , (3) 

n 

H ph = hw^blbn , (4) 
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Hint = ~ghuj ^ a n a n {b n + b n ) , 

■n 



(5) 



in which a' n creates an electron in the rigid-lattice Wan- 
nier state at site n, and b n creates a quantum hui of vi- 
brational energy in the Einstein oscillator at site n; J is 
the electronic transfer integral between nearest-neighbor 
sites, and g is the local electron-phonon coupling con- 
stant. 



II. CHARACTERISTIC ENERGIES 

The global ground state energy Eq is composed of con- 
tributions from the three principal components of the to- 
tal Hamiltonian: 

E = (ip\H\iP) = E km + E ph + E mt , (6) 

in which Ek in is the kinetic energy, E v ^ is the phonon 
energy, and Ei nt is the electron-phonon interaction en- 
ergy. 

Our numerical results for these quantities are summa- 
rized in Figure [l] as functions of the electron-phonon cou- 
pling strength g for several values of J/fiui. Each curve 
in Figure [l] contains 80-200 data points and each figure 
panel collectively represents nearly 1200 distinct polaron 
ground states. No smoothing has been performed; each 
"curve" is a polygonal arc connecting computed energies. 

Over most of the parameter space we have investigated, 
computational errors are smaller than can be meaning- 
fully conveyed with any resolvable symbol; the principal 
exception is at weak coupling and smaller J values where 
decreasing sensitivity of the ground state energy to cer- 
tain details of polaron structure eventually hampers con- 
vergence. 

With increasing adiabaticity, here beginning at 
J/hui « 7, the ability of the variational method to repre- 
sent the complexity of polaron structure in the immediate 
vicinity of the self-trapping transition eventually is over- 
taxed, and discontinuities appear in estimated quantities 
such as the kinetic energy. Although the value of the ki- 
netic energy in the immediate vicinity of such anomalies 
is necessarily distorted, the location of the discontinu- 
ities continues to provide reasonable estimates for the 
location of the self-trapping transition, and outside of a 
narrow region, quantitative accuracy remains good Jl8| . 
For such reasons, in the following we retain data for the 
cases J/hto = 7 and 9. 




FIG. 1. The ground state energy Eq (a) and its several components E^in (b), E p h (c), and Ei nt (d) as functions of the 
electron-phonon coupling for J/%u> = 0.25, 0.5, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0 and 9.0 (solid curves, top to bottom in (a), 
left to right in (b), (c), and (d)). The dashed curve in each figure panel is the exact J — > limit appropriate to each case. 
The chain-dotted curves illustrate Eqs. ((?]) - ( fll"| ) for J /Tiu = 9, and ( |l3| ) for J/hui = 1/4. Circles (o) indicate the "break" 
in each energy curve associated with the self-trapping transition. Bullets (•) in (a) and (b) indicate the J = termini of the 
self-trapping lines appropriate to _Eo and Ekm as determined from Eqs. (|ll|) and ( jl3| ) , respectively. The arrow in (b) indicates 
the initial slope of the E^ in self-trapping line as determined from Eq. ([L3p7 



Although the aggregate (Eq) of these several contri- 
butions is subject to variational constraint, the values 
of the separate contributions (Ekin, E p h, E p h) are not; 
in principle, the latter are vulnerable to distortions that 
may misrepresent some aspects of ground state structure 
while still yielding favorable results for the ground state 
energy. It is for such reasons that in the course of devel- 
oping our method and in using it to obtain new results, 
we pay close attention to the detail of overall polaron 
structure and compare multiple quantities with known 
results available from independent approaches. 
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In the appropriate regimes, the energy contributions 
shown in Figure [l] are in excellent agreement with 
the weak-coupling perturbation theory results (Rayleigh- 
Schrodinger, H = H krn + H ph , H' = H mt ) @-§| 



-23 - g 2 



Tilo 



E 



wc 



-23 + g 



E 



wc 

ph 



,/l + 4 3/hu ' 
, 23 
(1 + 4J/M 3/2 ' 
huj 23 



E 



wc 



y/1 + U/HUJ (1 + 4J/M3/ 2 

2%uj 
^/1 + 43/huj ' 



(7) 
(8) 
(9) 
(10) 



and with the strong-coupling perturbation theory results 
(Rayleigh-Schrodinger following the Lang-Firsov trans- 
formation, H = H ph + H int , H' = H km ) 
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where 7 is the Eulcr constant and Ei(y) is the exponential 
integral. One sample curve representing each of Eqs. (0) 
- (|l3]) has been included in Figure [l] to illustrate this 
agreement. More detailed discussion and comparisons 
with other theories have been given elsewhere [^8| . 

The kinetic energy is particularly important as a 
ground state property that is intimately connected with 
electron mobility. Characteristically, the kinetic energy 
is a weak function of the electron-phonon coupling below 
the self-trapping transition, and this dependence grows 
increasingly weak with increasing adiabaticity. Owing to 
the minimal involvement of phonons in the polaron in this 
regime, the quasi-particle can be fairly characterized as a 
quasi- free electron with a slightly reduced bandwidth. At 
very strong coupling, the kinetic energy decays to zero, 
suggesting that the dressed electron becomes essentially 
immobile relative to the quasi- free electron; significant, 
however, is the fact that this decay is not ultimately expo- 
nential in the coupling constant as is commonly assumed 
on the basis of the small polaron approximation (first 
order SCPT), but a much weaker inverse power as sug- 
gested by perturbative corrections (second-order SCPT) 
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A comparison of our numerical kinetic energies with 
weak (||) and strong ([l5| ) coupling perturbation theories 
is shown in Figure ||. Each of these comparisons exhibits 
sharp deviations in the vicinity of the self-trapping tran- 
sition because neither WCPT nor SCPT undergoes the 
transition while the variational results do. 




FIG. 2. Comparison of Global-Local kinetic energy, Ejf in , 
with weak-coupling perturbation theory, according to 

(^) (upper panel), and with the Lang-Firsov approximation, 
according to (|l^) (lower panel), as functions of the elec- 
tron-phonon coupling. 

Apart from data affected by deteriorating precision, 
the agreement between our variational and perturba- 
tive kinetic energy in the weak-coupling regime is very 
good up to g values close to the self-trapping transition. 
For this weak coupling case, perturbation theory char- 
acteristically overestimates the kinetic energy value, un- 
rferxepresenting the integration of phonons into polaron 
structure, causing the electron to appear more "free" 
than it actually is. This is to be expected, of course, 
since the inherent anharmonicity of the large polaron is 
induced by the electron-phonon interaction even at weak 
coupling, but is only partially captured by the low orders 
of weak-coupling perturbation theory. 

For strong coupling, our variational calculation ap- 
proaches the second-order SCPT result above the self- 
trapping transition; however, deviations persist signifi- 
cantly into the strong coupling regime, with larger J's 
converging to the SCPT result more slowly than smaller 
J's. Strong coupling perturbation theory systematically 
underestimates the kinetic energy, cwerrepresenting the 
integration of phonons into polaron structure, causing 
the electron to appear more "trapped" than it actually 
is. This, too, is to be expected, since the zeroth-order 
state of strong coupling perturbation theory is the ex- 
treme limit of a completely localized excitation; the finite 



G 



real-space spread of the true self-consistent polaron state 
is only partially captured at low orders. 




FIG. 3. Comparison of the kinetic energy Ekin as deter- 
mined by the Global-Local method (solid line) with weak cou- 
pling perturbation theory (dashed line), strong coupling per- 
turbation theory (dotted line), quantum Monte Carlo simula- 
tion (scatter plot), and density matrix renormalization group 
(diamonds) for J/Tilu = 1. The (+) symbol indicates the loca- 
tion of the self-trapping transition as determined here. QMC 
data kindly provided by P. E. Kornilovitch Q . DMRG data 
kindly provided by E. Jeckelmann . 

While this mutual consistency among limiting results 
is satisfying, the quantitative accuracy of both weak- 
and strong-coupling perturbation theory is superceded by 
that of the Global-Local method at intermediate coupling 
where the self-trapping transition is found. Figure || il- 
lustrates this point by comparing WCPT and SCPT with 
the Global-Local variational result as well as with results 
from quantum Monte Carlo simulation and the density 
matrix renormalization group approach. In the follow- 
ing, we rely exclusively upon our variational method for 
locating the transition. 



III. SELF-TRAPPING LINE 



The self-trapping transition is the more-or-less rapid 
change in polaron structure from that typical of large 
polarons (below the transition) to that typical of small 
polarons (above) as g or J are varied, and is typically 
evident in features that grow increasingly sharp with in- 
creasing adiabaticity. The self-trapping transition is cer- 
tainly the most exotic feature of the polaron phase di- 
agram, and is intimately involved with, if not always 
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ultimately responsible for, many of the difficulties en- 
countered in polaron theory. The physical transition is 
smooth at finite g and J/Tiu) |57|; however, as exempli- 
fied in the previous section, it is common for approximate 
descriptions of the phenomenon cither to miss the tran- 
sition completely, or to break down in some respect in its 
vicinity. 

In view of the smoothness of the physical self-trapping 
transition, however, there is no reason to expect the tran- 
sition to be manifested in exactly the same way in distinct 
physical quantities; thus, we expect an intrinsic ambigu- 
ity in the precise location of the the self-trapping transi- 
tion that is diminished only by the progressive sharpening 
of the underlying physical phenomenon. 

Since each of the energies displayed in Figure exhibits 
a "knee" or "break" that clearly separates distinct weak- 
and strong-coupling trends, any one of Eq 7 Ekin, E p h, or 
Ei n t could be used to locate the self-trapping effect. As 
objective criteria for locating this transition, we associate 
the transition with the particular coupling strength in 
the intermediate regime where each energy changes most 
rapidly with respect to g at fixed J/Tuo; these points are 
identified by zeros in the second derivatives with respect 
to g of Ekin, E p h, and Ei nt , and in the third derivative 
of Eq as collected in Figure 0; other rapidity criteria will 
produce slightly different traces. 
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FIG. 4. Polaron phase diagram, showing the location of 
the self trapping transition. Discrete symbols correspond to 
our numerical determinations by applying rapidity criteria to 
assorted Global-Local energy band data; we have included 
results from a similar analysis of the polaron effective mass in 
Ref. . The dashed line corresponds to the empirical curve 
gsr = 1 + \l J/hu; no fit has been performed to obtain this 



These criteria amplify numerical errors and there- 
fore are quite demanding of numerical precision. Over 
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the interval where the transition was sought, we used 
a sampling interval Ag sufficiently small that directly- 
computed g- values alone would suffice to convey the mes- 
sage of Figure || However, as our convergence crite- 
rion, wc required sufficient smoothness in the computed 
derivatives that the determining zeros could be interpo- 
lated with absolute numerical error considerably smaller 
than the sampling interval in most cases. While this de- 
gree of precision is not really necessary to convincingly 
locate the self-trapping line, we applied to our study 
the further "check" that the collection of numerically- 
determined self-trapping points should describe sensible 
curves through the actual family of energies as presented 
in Figure [l]. This proves to be a stringent error criterion, 
since the numerical errors already amplified by repeated 
differentiation with respect to g are further amplified by 
the characteristic steepness of Ekin, E p h, and Ei nt in 
the vicinity of the transition. It is especially significant, 
therefore, that the numerically-determined self-trapping 
loci obtained from Eki n not only describe a "sensible" 
curve, but appear to connect smoothly with the ascer- 
tainable trend as J/Tioj — > (See Figure |l|b). 

Such considerations permitted a clear resolution of the 
self-trapping transition for most hopping integral values 
used in this work, with a characteristic deterioration of 
precision at small J/Hlu and g here affecting primarily 
J/fruj < 1 and g < 1.5. For this reason, no self-trapping 
loci are reported for E p h and Ei nt at J/hcu = 1/4. 

Though precision deteriorates at small J/Huj, the dis- 
persion among estimated self-trapping locations that is 
evident in Figure [| at smaller J/huj values is not the re- 
sult of numerical errors, but reflects the intrinsic ambigu- 
ity in the assignment of a sharp location to a smooth tran- 
sition. Using perturbation theory (See (|i"T| ) and ([l3|)), 
for example, one can show that as J /hui — > 0, the trend 
in the kinetic energy criterion leads to g = ^= while 
the trend in the ground state energy criterion leads to 
g = *J\- This is consistent with the relative trends 
in the self-trapping estimates derived from Ej~i n and Eq 
data, and is at least roughly consistent with the absolute 
trends. This intrinsic non-coincidence implies that there 
is no single line consistent with all of the criteria that 
might legitimately be considered to locate the transition. 
Thus, the function ([l]) appears to describe the common 
or criterion-independent trend line about which distinct 
locators are tightly clustered with a spread that narrows 
steadily in both relative and absolute terms with increas- 
ing adiabaticity. 

There is reason to expect (|j) to continue to accurately 
locate the self-trapping transition beyond the investigated 
regime, since it is known from strong-coupling theory 
that the transition is associated with a particular value 
of the effective coupling parameter A = g 2 TiLo/2J in the 
adiabatic limit [01 



KsT = ~2J~ Xc 2 (16) 
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which is consistent with (|lj). 

Our conclusion that the relation (Q) accurately and 
robustly describes the location of the self-trapping tran- 
sition in the Holstein model is borne out as well in the 
behavior of other physical quantities not discussed here. 
Particularly important among these is the polaron effec- 
tive mass, which has been discussed in depth in Ref. Jl9| , 
whose consistency with our present results is illustrated 
in Figure [|. The effective mass is of great importance 
because it is a widely-recognized hallmark of the self- 
trapping transition and because it demonstrates that the 
scope of the self-trapping transition extends beyond the 
ground state. Other measures at the heart of the po- 
laron problem include the electron-phonon correlations 
that define the spatial structure of the polaron; presum- 
ably the self-trapping transition should involve charac- 
teristic rapid changes in spatial structure (e.g., "localiza- 
tion" upon self-trapping) that can be analyzed in order to 
assess the location of the transition. An extensive anal- 
ysis of such correlations to be presented elsewhere [|]]7]] 
fully supports the conclusions of this paper in quantita- 
tive detail. 



IV. ON DIMENSIONALITY AND 
ADIABATICITY 

The results of this paper apply to polarons in one space 
dimension. Extrapolating from these results to two and 
three dimensions is necessarily speculative; however, sev- 
eral pertinent observations can be made. 

The results that have long characterized commonly- 
held expectations for the dimensionality dependence of 
polaron structure are due to behavior ascertainable in the 
adiabatic approximation |3g)-|4C|]. Without recounting 
well-known arguments in detail, we note that these ex- 
pectations relate to the existence, possible co-existence, 
and relative stability of infinite-radius and finite-radius 
states; the former are often referred to as "free" states 
and the latter as "self-trapped" states. "Free" states are 
inessentially distinct from weakly-scattered free electron 
states while "self-trapped" states differ markedly from 
free-electron states in many respects. 

In two and three dimensions, the minimum energy 
states in the adiabatic approximation are found to be 
"free" states throughout the weak-coupling regime up to 
a discrete (structure-dependent) coupling threshold be- 
yond which "self-trapped" states have the minimum en- 
ergy. This transition phenomenon is what is meant by the 
term "self-trapping transition" in the adiabatic approxi- 
mation. Accordingly, there is no occasion to distinguish 
large from small polarons in two and three dimensions 
since the "free" states below the transition are of infinite 
radius and distinct from large polarons, and the "self- 
trapped" states above the transition are always inter- 
pretable as small polarons. This set of circumstances in 
two and three dimensions is reflected in the catch phrase 
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"all polarons are small" , since in this view large polarons 
in the adiabatic sense are never characteristic of the po- 
laron ground state. 

In one dimension, on the other hand, "free" states are 
not found at all in the adiabatic approximation except in 
the limit of vanishing electron-phonon coupling; instead, 
finite-radius (i.e. "self-trapped") states are found at all 
finite coupling strengths, leading to the commonly en- 
countered view that there is no self-trapping transition 
in one dimension. That these universally "self-trapped" 
states in one dimension might be distinguishable as large 
polarons or small polarons is of little consequence in this 
view, and one is led to consider the notion of a resolv- 
able transition between distinct large and small polaron 
structures as inconsequential as well. The issue cuts 
deeper, however, in that in the adiabatic approximation 
no such characteristic transition from large-to-small po- 
laron structure is found. 

That we here find clear and essential transition behav- 
ior between large and small polaron states in the one- 
dimensional case stands in stark contrast to the conven- 
tional adiabatic perspective; there is no casual reinter- 
pretation of one or the other set of results that relieves 
the contrast between such distinct alternatives. Since our 
one-dimensional results are quantitatively supported by 
independent h igh-qu ality methods (including cluster di- 
agonalization [OJ^-p4] , density matrix renormalization 
group [§5],|6) I an d quantum Monte Carlo p8|-|33f ) we are 
led to conclude that the problem to be resolved lies not 
with the data or methodology upon which we base our 
analysis, but with the adiabatic approximation itself. In- 
deed, we need not speculate on this point, since elabo- 
rations of adiabatic theory incorporating non- adiabatic 
corrections [^,^1 42 support our overall conclusion that 
the adiabatic approximation as it is widely regarded fails 
to embrace non-adiabatic characteristics that are essen- 
tial to the proper description of polaron states in the 
weak coupling regime, and therefore fails as well to prop- 
erly describe the self-trapping transition itself. This con- 
clusion is deeply rooted in the quantum-mechanical na- 
ture of electron-phonon correlations and is not unique 
to one dimension, suggesting that the adiabatic picture 
of self-trapping in two and three dimensions may mis- 
represent the true nature of polaron structure and self- 
trapping as well. We need not speculate too broadly on 
this point either, since an increasing body of results (es- 
pecially quantum Monte Carlo) consistently point to the 
absence of any dramatic qualitative differences between 
the occurrence of self-trapping in one dimension and in 
higher dimensions. 

The consistent conclusion from all these studies is that 
in any dimension self-trapping can be understood as a 
transition between large and small polarons, and that the 
large polaron state is non-trivially distinct from both the 
free electron state and the large polaron state commonly 
found in the adiabatic approximation. 

It may be suggested that although adiabatic correc- 
tions might be expected to be significant at sufficiently 
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small J/huj, the conventional adiabatic perspective on 
self-trapping ought to emerge as the adiabatic limit 
(J/hu> —> oo) is properly taken. As discussed at Jl6| ) 
ff., however, many arguments going beyond the adiabatic 
approximation point to a critical self-trapping transition 
in the vicinity of A c [l] = 1/2 in the adiabatic limit in one 
dimension, consistent with our findings in this paper and 
contrary to many widely-held expectations; i.e., our no- 
tion of a resolvable transition from large to small polaron 
structure persists into the adiabatic limit. 

Less firm information is available regarding the oc- 
currence of critical self-trapping in higher dimensions; 
however, in dimensions [D] dynamical mean field the- 
ory p3| , p| , for example, offers the estimates A c [l] = 1/2, 
A c [2] w 0.8, and A c [oo] « 0.844, while our own estimates 
based on scaling arguments [|7p^] suggest A c [l] = 1/2, 
A c [2] w 0.8536, A c [3] « 0.9082, and A c [oo] = 1. Both sets 
of estimates support the notion that at least the existence 
of the self-trapping transition and the scaling properties 
related to it are not seriously affected by changes in di- 
mensionality. 

Thus, adiabatic theory appears to be exceptional in 
suggesting a sharp distinction between polaron properties 
in one vs. higher dimensions. 

V. CONCLUSION 

In this paper, we have presented a large volume of 
data, comprehensive in scope and highly accurate, for 
the ground state energy Eq, kinetic energy Eki n , phonon 
energy E p h, and electron-phonon interaction energy Ei nt 
for the Holstein molecular crystal model of a single elec- 
tron in one dimension. 

We have broadly demonstrated the agreement between 
these results and the appropriate perturbation theories 
at weak and strong coupling, and have demonstrated the 
breakdown of both weak- and strong-coupling perturba- 
tion theory in the intermediate regime. For particular 
cases we have demonstrated the detailed mutual con- 
sistency of our results with those of the density matrix 
renormalization group and with those of quantum Monte 
Carlo simulation from weak coupling to strong coupling. 

Thus amply confirmed, we have analyzed the depen- 
dence of each of the principal energies upon the electron- 
phonon coupling constant to determine the self-trapping 
transition, defined relative to each measured property as 
the point at which that property experiences its most 
rapid change with respect to g. The data thus collected 
cluster in a clear way, suggesting the definition of a sin- 
gle self-trapping line cjst = 1 + \J J/^lu- This line joins 
with the adiabatic critical point in the adiabatic limit, 
and separates large polaron structure from small polaron 
structure at all J/huj. 

Whether one chooses to consider the specific data ap- 
propriate to one physical property or the gsr character- 
istic of all, the self-trapping points here determined are 
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consistent with features found in the particular kinetic 
energies, ground state energies, correlation functions, and 
effective masses determined independently by other high- 
quality methods. 

These findings support the notion of the self-trapping 
transition as a smooth phenomenon at finite J/huj and 
g, sharpening steadily as one approaches the adiabatic 
limit. As such, the self-trapping line associated with this 
transition is not unique, but is accurately and robustly 
located by gsr- 
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